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We study experimentally the propagation dynamics and interaction of a pair of mutually incoherent nematicons: spatial optical 
solitons in nematic liquid crystals. In contrast to earlier studies, we consider a bias-free liquid-crystal cell and compare the soliton 
interaction in copropagating and counterpropagating geometries. We analyze the dependence of nematicon interaction on input 
power and observe a direct manifestation of a long-range nonlocal nonlinearity. Attraction of counterpropagating solitons requires 
higher powers and longer relaxation times than that of copropagating nematicons due to losses-induced power asymmetry of 
counterpropagating nematicons. 
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1 Introduction 

Solitons have been observed in diverse fields of nonlinear physics, 
and they share common fundamental properties originating 
from the interplay between nonlinear self-action of wave packets 
and their natural tendency to spread as they propagate. Spatial 
optical solitons [1], i.e. nonspreading self- localized beams with 
the width unchanged during propagation, form due to a balance 
between linear diffraction and self-focusing in a nonlinear opti- 
cal medium. These solitons have been investigated extensively 
in several nonlinear media, both in one- and two-dimensional 
geometries; they have significant potentials in signal processing, 
switching and readdressing in the future generation of all-optical 
circuits. In this context, a giant optical nonlinearity arising from 
molecular reorientation in nematic liquid crystals (NLCs) has 
attracted significant attention [2-4]. Both experimental [5,6] 
and theoretical [7] results have been demonstrated for spatial 
optical solitons in nematics, also called nematicons [8]. 

Nematic liquid crystals consist of elongated molecules aligned 
along a given direction (known as the molecular director) ow- 
ing to both anchoring at the boundaries and intermolecular 
forces [2-4]. The resultant medium exhibits a positive uniaxial 
anisotropy and birefringence, with ordinary and extraordinary 
refractive indices, nu and n±, defined for polarizations parallel 
and orthogonal to the director. The reorient at ional nonlinearity 
allows generating nematicons at relatively low optical powers, in 
the milliWatt region or below, for the study of the fundamental 
aspects and applications of light interaction with self- assembling 
nonlinear soft matter. 

While initial studies of solitons in nematic liquid crystals con- 
sidered the propagation of single nematicons [5,9-14], the in- 
teraction of two copropagating nematicons of equal [15-19] or 
different wavelengths [20,21] have been addressed more recently. 




Figure 1: Top view of the NLC planar cell for the study of (a) copropagat- 
ing and (b) counterpropagating solitary beams. The ellipses indicate the 
orientation of the molecular director in the plane (y, z). (c) Experimental 
setups for the copropagating (solid lines) and counterpropagating (dashed 
lines) cases: Ls - cw-laser (A = 532nm), BS - beam splitters, M - mirrors, 
A/2 - half wavelength plates, MO - microscope objectives, NLC - cell, F 
- filter, MS - microscope, CCD - camera. 

The attraction and fusion of nematicons in bulk undoped NLCs 
in planar voltage-biased liquid crystal cell were reported by Pec- 
cianti et al [15,22]. 

Most of the studies considered nematicons propagating in the 
same direction. At the same time, the interaction of counter- 
propagating solitons in, e.g., photorefractive crystals [23-26], 
is known to posses rich physical behavior, including convective 
dynamical instabilities [27-30]. The latter were observed in ex- 
periment in Ref. [31] and stabilization against these instabilities 
was achieved in transversely periodic nonlinear medium [32]. In 
contrast, there is little known about counterpropagating nemati- 
cons, although an experimental approach was developed, based 
on the mutual deflection of two counterpropagating beams, 
which permitted to estimate the strength of thermal focusing in 
thick dye doped nematic-liquid-crystal samples excited by nar- 
row laser beams [33,34]. In addition, recent numerical studies 
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Figure 2: Experimentally recorded images of light scattered from two CO 
nematicons at various input powers P in each beam. 

pointed out the possible occurrence of instabilities for counter- 
propagating beams in liquid crystals [35,36]. 

In earlier experimental studies with nematicons the liquid- 
crystal cell was biased by an external electric field, aimed at con- 
trolling the orientation of the NLC molecules at rest and, there- 
fore, their nonlinearity as well as their nonlocality [37] , thus also 
affecting nematicon interactions. However, nematicons can exist 
in unbiased NLC cells provided the input beam is extraordinary 
polarized, e.g. with director and electric field coplanar with the 
plane (x,z) parallel to the cell interfaces [13, 14, 17, 18,20,38]. 
Because of the optical anisotropy of liquid crystalline molecules 
and birefringence, light beams propagating in NLC walk-off the 
direction of their wave vector. The walk-off can be adjusted by 
acting on the optic axis, i.e. by reorienting the molecular direc- 
tor [13,39]. 

In this paper we experimentally address the propagation and 
interaction of copropagating (CO) and counterpropagating (CP) 
nematicons in a bias-free NLC cell. We investigate in detail both 
power- and time-dependent interactions of two identical solitons. 
In particular, we demonstrate that the long-range attraction and 
ability of nonlocal solitons to form bound states are strongly 
affected by the geometry of the interaction. 



2 Experimental results 

We employ the cell geometry sketched in Fig. E{a), with two 
parallel polycarbonate plates separated by a 100/xra gap. The 
cell contains the 6CHBT liquid crystal [40, 41] which has neg- 
ligible absorption and high nonlinearity with refractive indices 
n e = 1.6718 and n — 1.5225 at room temperature. The polycar- 
bonate plates are rubbed in the plane (x,z) at the angle of 
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Figure 3: Experimental results of the attraction between two CP solitons 
at different input powers, (a) Measured distance r between solitons (at 
propagation distance z — 0.55 mm) and corresponding (b) normalized 
intensity profiles for various input powers in each beam. 

45° with respect to the z axis. Such prepared surfaces entail 
molecular orientation in a bulk analogous to that provided by 
pre- alignment via an external biasing field when the director is 
aligned along z. However, in contrast to the latter case where 
the director can rotate in a plane orthogonal to the plates, here 
the optical axis remains parallel to them. 

The experimental setup is shown in Fig. \Hc). Two identical 
beams generating spatial solitons are prepared using a system 
of mirrors and beam splitters. For CO solitons, two parallel in- 
put beams of extraordinary polarization (E-fleld along the x- 
axis) are formed using a standard Mach-Zehnder arrangement 
(beamsplitter BS1 and BS2, mirrors Ml and M2). In the case 
of CP solitons, the mirror Ml is removed and the beam trans- 
mitted by BS3 follows the route indicated by the dashed line. 
In addition, in the case of CO beams, the mirror Ml mounted 
on a piezoelectric transducer is made to vibrate (at frequency 
1kHz) to induce a frequency shift on the reflected beam. Since 
the nonlinear response of NLC is slower than the mutual phase 
changes of two beams ~ 10 -3 s, such frequency detuning makes 
the interaction of beams effectively incoherent. With the half- 
wave plate we control the polarization state of the two identical 
and parallel Gaussian beams focused into the cell by a 10 x ob- 
jective. The waist of each beam is wo = 2 /am and the distance 
between them is 25 fim. The beam evolution inside the cell is 
monitored by collecting the light scattered above the cell with 
a 5x microscope-objective and a high resolution CCD camera. 
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The total beam power in all experiments was kept low (below 
15 mW) in order to avoid thermooptical effects arising from 
light absorption and heating. 

First study the dynamics of the formation and interaction of 
two CO parallel beams in NLC as a function of the excitation 
power, see Fig. [2] As discussed earlier [15], due to of the long- 
range character of nonlocality, the nematicons attract each other 
progressively more for increasing input power. In most cases 
this attraction independent of the relative phase of the solitons, 
even when they are coherent to each other. A sequence of im- 
ages in Fig. [2] shows the stationary trajectories of two initially 
parallel beams for different input powers. For low input power 
(0.5 mW), the self- focusing is too weak to overcome diffraction, 
so the beams keep spreading as they propagate. By increasing 
the power to 2 mW, we achieve a stable propagation of soli- 
tons and their weak attraction. Due to high birefringence of the 
liquid crystal, both nematicons propagate at an angle with re- 
spect to the initial wave vector (directed along the z-axis). The 
measured walk-off of the Poynting vector is ~ 4.3°. For higher 
input power the attraction is so strong that it induces one (at 
P > 3 mW) or multiple (P = 13.5 mW) intersections of the 
soliton trajectories. 

Next we investigate the power-dependent dynamics of the CP 
solitons in the same cell. In this case the beams are launched 
from the opposite sides of the cell by focusing them with 10 x 
micro-objectives, see Figs. [I{b,c). F° r better comparison with 
the CO case discussed above, we align the initial trajectories 
(without walk-off and without interaction) to be parallel to the 
optical axis z and spatially separated by 25 /mi. As in the case 
of CO nematicons, the waist of each CP beam is estimated to 
be wq = 2 /im. 



Experimental results for CP beams are presented in Fig. [3] As 
above, increasing the input power (P >2 mW) leads to self- 
focusing and the formation of two CP nematicons. However, by 
further increasing the power we observe a clear difference in the 
trajectories of the CP nematicons compared to the CO ones [cf. 
Fig. [2] and Fig. [3]. In particular, we find that the CP trajecto- 
ries remain almost parallel up to a power of 5 mW while the 
strong attraction in the CO case is already observed for about 
2.5 mW. For higher powers, the trajectories approach each other 
until they eventually merge into a single soliton at 13.5 mW. In 
this regime both solitons propagate along the same trajectory 
forming the so-called vector soliton. This behavior is illustrated 
in Fig. [3jb) where we plot the separation between trajectories 
of CP nematicons versus excitation power. It is clear that the 
trajectories start approaching each other only after the input 
power exceeds 5 mW. For higher powers, up to 18 mW, we do 
not observe changes in the soliton positions, neither associated 
with crossing, as in the CO case, nor because of spatiotemporal 
instabilities. Figure [3jc) shows the intensity profiles for different 
input beam powers at the propagation distance of z = 0.55 mm. 

In order to understand the origin of the different dynamics of 
CO and CP nematicons, we have to consider the role of scatter- 
ing losses, which progressively reduce the beam power and the 
size of the nonlinear effect. While the two CO beams evolve for- 
ward with equal individual powers at each z, the CP beams 
interact with unequal intensities and sizes, i.e. with unequal 
strengths in the transverse force pulling the two nematicons to- 
wards each other. The two initially parallel CP solitons, say A 
and B, launched at z = and z = L, respectively, tend to at- 
tract and shift the center of mass of the resulting dipole towards 
the more intense one, i.e. towards A for z < L/2 and towards 
B for z > L/2. This in- homogeneous force distribution results 
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Figure 6: (a) The bound CO nematicon pair at t = for the power of 
13.5 mW and the initial separation of 25 /im (the same as in bottom frame 
of Fig. [2}. (b, c) Superimposed intensity profiles of a single beam in the 
pair (for t — 0) and for t = 5 s after the second beam has been blocked. 
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Figure 7: The same as in Fig.[6]for CP nematicons from the bottom frame 
of Fig. [3] The white line in (b) and (c) shows the position on the intensity 
maximum of the coupled pair in (a). 

in a larger power requirement for the completion of the merging 
process as compared to the CO case. 

The asymmetry between CP nematicons has a profound effect 
on the temporal dynamics of interaction inducing significant 
differences with the CO case. Experimental sequence of images 
from Video 1, demonstrating the temporal evolution of two ini- 
tially parallel CO beams, is presented in Fig. 2] These images 
clearly reflect the inertia of the nonlinearity. After the light is 
switched on at t = s, the beams propagate linearly and diffract. 
It takes some time (~ 0.3 s) for the solitons to form. Only after 
additional time the soliton interaction becomes visible, leading 
to multiple intersections of their trajectories. In a sharp con- 
trast, the temporal evolution of the CP beams is much slower 
as it requires 17 s for the solitons to merge and form a bound 
state, see Fig. 2] and Video 2. 

We stress here that the differences between CO and CP nemati- 
cons are solely due to their interaction as all individual nemati- 
cons are equivalent in our experiments. To demonstrate this fact 
we study the relaxation dynamics of each beam after its neighbor 
is switched off at t = 0. The results for CO (CP) nematicons 
are presented in Fig. [6] (Fig. [7]). To highlight the differences, 
in Figs. EJb,c) we superimpose two images of the same beam, 
with and without interaction, at t = and t = 5 s respectively. 
Similarly, the two "free" CP nematicons in Figs.[7Jb,c) are com- 
pared to their common waveguide (white line) . In both cases the 



figures demonstrate how the remaining soliton reconstructs its 
original individual trajectory after about 5 s. Importantly, this 
relaxation time is the same for both CO and CP cases because 
it is spent by a single nematicon to erase the changes in director 
orientation induced by both interacting beams. Thus, an order 
of magnitude difference in times necessary to build bound states 
in Figs. H] and [5] should be considered as an essential signature 
of symmetric vs. asymmetric interactions of nematicons. 

3 Conclusions 

We have investigated experimentally the generation, temporal 
dynamics, and power-dependence of interacting copropagating 
and counterpropagating spatial solitons in unbiased nematic liq- 
uid crystals outlining differences between these two cases. In 
particular, we have observed that the attraction of copropagat- 
ing nematicons occurs at lower powers than for counterpropa- 
gating solitons, with an additional significant difference in the 
time response between the two cases. 

This work was supported by the Australian Research Council. 
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